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1. Introduction 

In the interface region of coupled ferro- and antiferromagnetic systems there is a 
magnetic reordering known as the magnetic proximity effect (MPE). MPE has for a 
long time attracted the interest of researchers [1] since the constituents show a novel 
magnetic arrangement different from that of the bulk. The interest in such interfaces 
has revived lately with respect to the exchange bias effect [2] , which occurs when a thin 
ferromagnetic (FM) film is deposited on an antiferromagnetic (AFM) material. If the 
latter is an 'in-plane AFM', and the number of bonds between parallel and antiparallel 
spin pairs across the interface is the same, the interface is 'compensated'. In this case 
the AFM often assumes an almost orthogonal magnetization with respect to the FM 
magnetic direction, while the spins of the AFM interface layer reside in a 'spin-fiop- 
phase', in analogy to an AFM system in an external magnetic field [3]. 

In previous studies investigating FM-AFM interfaces the magnetization of each 
FM layer is usually considered to be coUinearly ordered and to rotate as a whole 
[4, 5]. Results concerning the spin reorientation transition (SRT) have been obtained for 
various magnetic systems but. to the best of our knowledge, only in presence of magnetic 
anisotropics [5]. It is important to stress that, although we consider in the present work 
anisotropic interactions too, a noncoUincar magnetization in coupled FM-AFM bilaycrs 
is caused mainly by isotropic exchange interactions. Whereas previous work describes 
a net magnetic reorientation of the total system, we consider the reorientation to take 
place in magnetic sublattices, leaving the net magnetic orientation of each layer almost 
unchanged. In the case of a compensated FM-AFM interface, the MPE extends to 
the FM layers close to the interface. Then, the magnetic structure of each FM layer 
is represented, in perfect analogy to the AFM layers, by two juxtaposed sublattices 
with different but uniform magnetization directions. Allowing a nonuniform intralayer 
magnetic structure in the FM subsystem leads to new features, which in turn are strongly 
dependent on the underlying lattice symmetry. 

Often for simplicity, a mean field theory (MFT) is used to describe the magnetic 
reorientation. In a recent MFT [6] study of a coupled FM-AFM system for both sc(OOl) 
and fcc(OOl) lattices, a variety of different magnetic configurations emerge, depending 
on the parameter values. Usually the subsystem with the larger ordering temperature 
induces a magnetic order in the other one (MPE). For coupled sc(OOl) systems, both 
FM and AFM films are perturbed from their collinear magnetic order and exhibit 
similar behavior. This symmetry is absent for fcc(OOl) films, which, under certain 
circumstances, may exhibit two different critical temperatures. An advantage of MFT 
is that many results can be derived analytically for simple bilayer systems. 

We stress that for thin magnetic films, collective magnetic excitations (spin waves) 
are particularly important. These are neglected completely by MFT. Such excitations 
are taken into account, for example, by many-body Green's function theory. There is a 
large amount of work applying this theory to thin ferromagnetic films, in particular in 
connection with the reorientation of the magnetization as a function of temperature and 
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film thickness. We mention only a few papers that cite further literature [7, 8, 9, 10]. 
Antiferromagnetic films have also been treated, e.g. in [11, 12]. Not as much work has 
been done in which Green's function theory treats the coupling of ferromagnetic layers 
to antiferromagnetic layers: in reference [13], a bilayer is investigated and reference 
[14] treats an extension to multilayers. In both cases, only a coUinear magnetization is 
considered. In reference [15], a ferromagnetic film is coupled to an antiferromagnetic 
layer; however, the orientation of the magnetization of the antiferromagnet is frozen. 
Other work considers an antiferromagnetic coupling between ferromagnetic layers 
[16, 17, 18]. 

The new feature of the present work is to allow an in-plane reorientation of the 
magnetizations of the ferro- and antiferromagnetic layers due to the interlayer coupling 
as in the MFT approach of [6] but using Green's function theory. We restrict ourselves 
to Heisenberg systems with spin S = 1/2 with an exchange anisotropy. We do not 
consider this an essential restriction, because we have shown in references [19, 20] for 
ferromagnetic layers that through an appropriate choice of anisotropy parameters the 
exchange- and single-ion anisotropies yield very similar results, and that an appropriate 
scaling leads to universal magnetization curves for different spin quantum numbers. 
In the present paper, we examine in detail the magnetic arrangement of the simplest 
system: a perfectly ordered bilayer consisting of a FM monolayer that is coupled to a 
AFM monolayer. Thus, we do not address the exchange bias effect directly, since this 
effect is most likely related to a certain amount of interface disorder [21]. 

The paper is organized as follows. In Section 2 we develop the Green's function 
theory and discuss the method of solution for the resulting equations. In Section 3 we 
present numerical results for bilayers with purely ferromagnetic or antiferromagnetic 
couphng as well as for the coupled FM-AFM bilayer. In particular the effect of the 
interface coupling Jint on the characteristics and magnitude of the MPE at zero and finite 
temperatures is investigated. The resulting magnetic arrangements for various kinds of 
bilayer systems and for their corresponding ordering temperatures are determined. In 
Section 4 we summarize the essential results and end with some remarks concerning 
further development. 

2. The Green's function formalism for coupling ferro- and 
antiferromagnetic layers 

The starting point is a XXZ-Heisenberg Hamiltonian consisting of an isotropic 
Heisenberg exchange interaction with strength J^j between nearest neighbour lattice 
sites, exchange (non-localized) anisotropies in the x- or z-directions having strengths 
Df^ and Df^ respectively, and an external magnetic field B = [B^ B'^) confined to 
the film plane, which is the a;2;-plane: 
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<ij> <ij> 

-Y.{B^S^ + B^Si). (1) 

k 

Here the notation 5*^^ = Sf ± iS^ is introduced, and < ij > indicates summation over 
nearest neighbours only, where i and j are lattice site indices. Because there is no field 
perpendicular to the film plane, only a reorientation of the magnetization in the 
x^-plane is allowed. For the FM-AFM bilayer we use D^j in the ferromagnetic layer and 
Dfj in the antiferromagnetic layer. 

In this paper we restrict ourselves to systems with spin quantum number S — 1/2 
and to a simple cubic (sc) lattice. For this case we need the commutator Green's 
functions 

Gt-{u;) = {{Sr;S-))^, (2) 

where a = (+, — , z) takes care of all directions in space. A generalization 
to spin quantum numbers S > 1/2 is straight- forward by introducing G"''"" = 

; with m + n < 2S + 1 (m > 0; n > 1; m,n integer) [7]. 

The equations of motion for the Green's functions in the energy representation are 

u;Gtr{u) = A^ + {{[Sr,nSi)). (3) 
with the inhomogeneities 

( 2(5f)5., \ 

^f-{[Sr,S.])= , (4) 

V -{Sf)k J 

where (...) = Tr(...e^^''^)/Tr(e^'^^) denotes the thermodynamic expectation value. 

In order to obtain a closed system of equations, the higher-order Green's functions 
on the right hand sides arc decoupled by a generalized Tyablikov- (RPA) decoupling [7] 

{{SrS^k-, Si)U 0, {St)Gij + {Sl)Glr. (5) 

After introducing two sublattices per layer, which is necessary when dealing with 
antiferromagnets, the resulting equations are Fourier transformed to momentum space 

by 

fllfl \ / AT / J (■mjn ' 

GZjn =|^EGrn(k)e^''(^--^-\ (6) 

where im, jn are lattice site indices on the sublattices m,n, and N is the number of 
lattice sites in the whole system. One obtains 
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T (5^ + E('5;)('^-p(o) + ^^p(o)))G^n 
p 

±(^^)E(-^-p(k)+L'^,(k))G'-, 

p 

^ p 

+ 9 ('S'm) E ■^"iplkjG'pn 
^ P 

+ k^-^ + E('^;)(>^™p(o) + i^^p(o)))g-; 
p 

-^('5;?.)E'^-p(k)G,T- (7) 
^ p 

For a square lattice with lattice constant qq — 1, one has four nearest-neighbour 
intmlayer couplings with sublattice indices n, m from the same layer 

<^mn(0) = Q'o Jmn ■> <^mn(k) = 7o(k) Jmn i 

(8) 

D^^m = go D^^l , D;^^(k) = 70 (k) D^^l , 

with the intralayer coordination number qo = 4 and the momentum-dependent Fourier 
factor 

7o(k) = 2(cosA;a; -|- cos kz) ■ (9) 

Correspondingly, for the nearest neighbour interlayer couplings, with m, n now being 
sublattice indices from different layers, one obtains 

JmniS^^ — Qint J'mt i <^mn(k) — 7int(k) Jint i 
^^rn(0) = ?int A^n^ , ^^^(k) = 7int(k) A^J^^- 

For sc stacking, the m^eWa?/er coordination number and the corresponding Fourier factor 
are given by 

9int = 7mt(k) = 1 , (11) 

which is assumed in the following calculations if not stated otherwise. 
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For fee or bee staeking one has for eomparison 

gint = 4 and 7i„t(k) = 4 eos(A;^/2) cos{k^/2). (12) 

The mean field approximation is obtained by negleeting the Fourier factors, i.e. 70 (k) = 
7int(k) = 0. By ehoosing the appropriate signs of the exchange interaction and the 
exchange anisotropy couphng constants, one can treat ferromagnetic, antiferromagnetic, 
and mixed systems with coupled FM and AFM layers. 

The general formalism is valid for any number of layers and sublattices. If Z is 
the total number of sublattices of the system, the dimension of the set of equations 
(7) is 3Z^. Because we restrict ourselves in the present paper to the investigation of 
the bilayer problem, there are four sublattices, and the system of equations (7) is of 
dimension 48 with a corresponding Green's function vector. Closer inspection reveals 
that the system of equations has the following substructure 

/ /rooo\\/Gi\ /Ai\ 



ujl - 



A2 
A3 

V A4 y 



(13) 



r G2 
r G3 

V V r yy V G4 y 

where the diagonal blocks T are identical 12 x 12 matrices, whose explicit form can be 
read off from equations (7). The sublattice Green's functions G„ {n — 1,2,3,4) are 
vectors of dimension 12 consisting of 4 subvectors, each of dimension 3: 

/ Gi„ \ 

G2n 



y 



Gan 

V Gin 

where the 3-component vectors are 

/ G+- \ 

m.n 



71 = 1,2,3,4, 



(14) 



G. 



\ G^ 



m 



1,2,3,4. 



(15) 



y 



All, — 



A — 



( 2(5^) \ 




■m,n — 1,2,3,4 . 



(16) 



The inhomogeneity vectors have the same structure: 

/ Ai„(5i„ \ 

A2n^2n 
Asnf^Sn 
V ^An^An / 

The big equation (13) of dimension 48 for the bilayer can therefore be replaced by 4 
smaller equations of dimension 12: 

(a;l-r)G„ = A„ for 7i = 1, 2, 3, 4 . (17) 

By invoking the spectral theorem, one obtains an expression for the diagonal 
correlations in configuration space C„ (independent of k) by integrating over the 
correlations in momentum space C„(k) corresponding to the Green's functions G„ [7]: 



CJk) 



(ikRsLA^ 



71 = 1,2,3,4, 



(18) 
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where R and L are matrices constructed from the right and left eigenvectors of the 
eigenvalues of the non-symmetric matrix F, and £ is a diagonal matrix with elements 
Sij/^e^"^^ — 1) obtained from eigenvalues uji {i = 1, 12) in the case when all of them 
are not zero. 

Unfortunately, equation (18) cannot be used directly because the 12 x 12 F-matrix 
turns out to have 4 zero eigenvalues. These have to be treated properly when applying 
the spectral theorem in order to calculate the correlations. In reference [7] we have 
shown that in this case the correlations C„(k) obey the equations 

(l-R°L°)C„(k) =RVL^A„, 71=1,2,3,4, (19) 

where 



R — (R'^j •^^) ) 




The subscripts and 1 label the matrices constructed from the right and left eigenvectors 
belonging to the zero and non-zero eigenvalues respectively, and is a diagonal 
8x8 matrix with elements 5ij/{e^'^^ — 1) obtained from the non-zero eigenvalues 
uJi {i = 1, ...,8). 

As written, equation (19) is formally correct but it still cannot be used directly for 
two reasons: (1) The matrix (1 — R'^L^) is a projection operator on the non-null space 
and thus has no inverse; i.e. we cannot isolate C„(k) unless (1 — R°L°) is independent 
of k. In the example of reference [22], this is the case and one take the projector outside 
the integral. This is not the case in the present example. (2) Even if (1 — R°L°) were 
independent of k, we still do not know how to treat the non-diagonal correlations Cmn 
of a multilayer problem. 

We have considered these problems in recent publications [22, 23], where we show 
that a solution is attained via the the singular value decomposition (see e.g. [24]) of the 
F-matrix 

F = UWV = uwv, (21) 

where U and V are orthogonal matrices, and W is a diagonal matrix with the singular 
values on the diagonal, and u and v are matrices obtained from U and V by omitting 
the colums and rows corresponding to singular values zero. Multiplying equation (19) 
from the left by v and inserting vv -|- VqVo = 1 gives 

v(l - R°L°)C„(k) = RV1L^(vv + vovo) A„ , (22) 

and use of vR° = 0, LVq = 0, and r = vR^ and 1 = L^v leads to the equations 

vC„(k) = r£ilvA„ , 71=1,2,3,4. (23) 

For more details of the formahsm consult reference [23] , where we present a systematic 
way of finding for each sublattice one k-independent v vector having a layer structure, 
i.e. v„ = (0, 0, .., 0). In this way the non-diagonal correlations disappear from 
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those rows in equation (23) corresponding to v„, and the k-integration can be performed: 
/ cikvC„(k)= V / cikC„(k) = vC„. In the present case v„ is given by 



Putting equation (24) in equation (23) yields 4 equations which contain the 8 
magnetization components imphcitly. The necessary additional 4 equations are obtained 
from the regularity conditions [7] 



which are obtained from the fact that the commutator Green's functions must be regular 
at the origin, see e.g. [25]. The Uq are constructed from the null-eigenvectors of the 
singular value decomposition of the matrix T. The resulting 8 integral equations are 
solved self-consistcntly by the curve-following method described in detail in the appendix 
of reference [26]. Note that the uq are determined numerically only up to an orthogonal 
transformation. To ensure proper behaviour as a function of k, iio must be calibrated 
at each k. A procedure for effecting this is presented in an appendix of reference [23] . 

3. Results 

In the following we present results for the bilayer ferromagnet, the bilayer 
antiferromagnet, and the coupled ferro- and antiferromagnetic bilayer. All calculations 
are for an in-plane orientation of the spins of both layers. In each case we compare 
the results of Green's function theory (GFT) with those of mean field theory (MFT) 
obtained by putting the momentum-dependent terms equal to zero. 

In order to see the effects of the interlayer coupling most clearly, we use different 
exchange interaction strengths for both layers: 

(a) FM-FM: Jifm = 100, Jsfm = 50, 

(b) AFM-AFM: Jiafm = -100, J2AFM = -50, 

(c) FM-AFM: Jfm = 100, Jafm = -50. 

Because of the Mermin- Wagner theorem [27], one needs anisotropics in the Green's 
function description: We assume — -1-1.0 for FM layers and — —1.0 for AFM 
layers. The magnitude of the anisotropies is appropriate for 3d transition metal systems. 
For a compensated interface, the magnetizations of the FM and AFM layers are almost 
orthogonal to each other even at T = because of the interface exchange interaction 
-^int- We choose the FM magnetization to be oriented in the 2;-direction and the AFM 
magnetization in the ,x-direction. Our particular choice of the anisotropies supports 
this arrangement not only at T = but also at finite temperatures. For other 
choices of anisotropies the magnetic arrangement could be different. The interlayer 
coupling is assumed to be positive for the ferromagnetic bilayer and negative for the 
antiferromagnetic bilayer. For the coupled FM-AFM system, both signs are used. 




(24) 




n = 1,2,3,4 



(25) 
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We consider three interlayer coupling constants with strength Jint = 30, 75, 160, 
respectively, one smaller than the weakest exchange interaction, one larger than the 
strongest exchange interaction, and one in between. 

3.1. The ferromagnetic and the antiferromagnetic hilayers 

Results for the FM and AFM bilayers are presented in this subsection. We do so in order 
to have a basis for discussing the differences to the coupled FM-AFM bilayer described 
later on. 

In figure 1(a) we show the sublattice magnetizations of the ferromagnetic bilayer 
as a function of the temperature for three interlayer couplings calculated with Green's 
function theory (GFT). The magnetization profiles are different for the two layers (the 
magnetization is larger for the layer with the larger exchange interaction) but end in a 
common Curie temperature, which increases with increasing strength of the interlayer 
coupling: Tcurie = 50.66, 55.24, 60.04. 

For the antiferromagnetic bilayer we use the same parameters as for the 
ferromagnetic bilayer except for a sign change. In figure 1(b) we show the sublattice 
magnetizations of the antiferromagnetic bilayer for two interlayer coupling strengths 
calculated with Green's function theory. To make the figures more transparent we leave 
out the result for the intermediate interlayer coupling strength. The corresponding 
magnetization curves lie in between those of the other couplings. One observes clearly 
the well known reduction of the magnetizations at low temperatures due to quantum 
fluctuations, which are missing in MFT, see figure 1(c). Since |Jiafm| > |-^2Afm| this 
reduction is larger for the first layer. With increasing temperature the magnetization 
curves of the two layers cross each other, a fact which was first observed by Diep [11], 
and finally end in a common Nccl temperature. A larger interlayer coupling leads to 
a larger suppression at low temperatures and to a larger Nccl temperature. Whereas 
with the present choice of parameters the magnetization profiles of the FM and AFM 
bilayers are rather different at low temperatures, the critical temperatures turn out to be 
identical: Tcurie — ^N^ei (cf figures 1(a) and (lb)), a fact that has already been discussed 
by Lines [28]. 

For comparison, we show in figure 1(c) the results of mean field theory (MFT) with 
the same parameters. The magnetization profiles as well as the critical temperatures 
are identical for the ferromagnetic and antiferromagnetic bilayers. As is well known, 
the Curie (Neel) temperatures (Tcurie(Neei) = 102.10, 107.25, 123.16 ) are much larger in 
MFT due to the missing magnon excitations, with the present choice of parameters by 
about a factor of 2. Note that in MFT the Curie temperature is not very sensitive to the 
anisotropics, as long as they are much smaller than the exchange interaction. In GFT, 
however, the sensitivity is very much greater because of the Mermin- Wagner theorem 
[27] (Tcurie(Neei) for D^(^) — > 0). One observes also that the effect of the interlayer 
coupling on the magnetization profiles is much stronger in MFT than in GFT. 
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Figure 1. (a) Green's function theory (GFT) for the ferromagnetic bilayer: The 

sublatticc magnetizations are displayed as a function of the temperature for different 
interlayer couplings Jint = 30 (dotted), 75 (dashed), 160 (solid). The exchange 
interaction and anisotropy constants are Jifm = 100, J2FM = 50, D^ym — 
1.0, D^^, = 1.0. 

(b) GFT for the antiferromagnctic bilayer: The sublattice magnetizations are 
displayed as a function of the temperature for two interlayer couplings Jint = 
—30 (dotted), —160 (solid). The exchange interaction and anisotropy constants are 

JlK¥M = —100, J2AFM = —50, D^p^YM = ^l-0j-O2AFM = —1-0. 

(c) Mean field theory (MFT) for the ferromagnetic and antifcrromagnetic bilayers 
with identical parameters: Ji(2)fm = ki(2)AFM|> -Di(2)fm = 1-Di(2)afm|> -'intFM = 

l-^intAFMl- 



3.2. The coupled ferro- antiferromagnetic bilayer 

This is the most interesting case. In the present study, we consider two in-plane 
magnetization components of each sublattice thus allowing noncoUinear magnetizations 
in both the FM and AFM layers. Our computer code, when speciahzed to a single 
magnetization direction, reproduces the results of reference [13]. Without interlayer 
coupling, the code also reproduces the results for both the monolayer ferromagnet and 
monolayer antiferromagnet simultaneously. The present choice of anisotropies supports 
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Figure 2. (a) Green's function theory(GFT): The sublattice magnetizations of the 
ferro- and antiferromagnetic sublattices are displayed as a function of the temperature 
for different interlayer couplings Jint = 30, 75, 160. The exchange interaction and 
anisotropy constants arc Jfm = 100, Jafm = —50, -D^pj^^ = —1.0, -Df-M = 1-0. 
(b) Mean field theory (MFT): with the same parameters. 



the orthogonal arrangement of the magnetizations of the FM and AFM layers favoured 
by the exchange interaction alone. The interlayer coupling destroys the perpendicular 
orientation of the ferromagnet (in z-direction) with respect to the antiferromagnet 
(in x-direction) even at temperature T = 0, as can be seen from figures 2 and 3. 
In figure 2(a) we show the sublattice magnetizations of a sc lattice calculated with 
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GFT for three interlayer coupling strengths. With a positive interlayer coupling all 
sublattice magnetizations develop a positive z-component, whereas the x-components of 
the two sublattice magnetizations in each layer are oriented oppositely. With increasing 
temperature, all x-components decrease until they vanish at a common temperature 
^Neeb slightly above the Neel temperature of the uncoupled AFM. For T > Tj^^gi all 
sublattice magnetizations point in the positive z-direction, the AFM layer assumes a 
ferromagnetic arrangement, and remains so until a common critical temperature Tq is 
reached, at which the magnetic order vanishes altogether. This is more clearly shown 
in figure 3 for the case of the strongest interlayer coupling of figure 2. Due to the 
strong interlayer coupling, the magnetizations of the layers are no longer coUinear, even 
at T = 0, and turn more and more into the ^-direction with increasing temperature 
until T^(5ci! while the magnitudes of the magnetization vectors shrink. Above T^^^j the 
sublattice magnetizations stay coUinear until they vanish at the critical temperature Tc- 



FM and AFM sublattice magnetizations 



z t 



FIVI 



AFM 





-10 



10 

temperature 



20 



Ti 



30 



Neel 



40^ 



50 



Figure 3. The reorientation of the sublattice magnetizations of the FM and AFM 
bilayer as a function of the temperature for the case of the strongest couphng 
( Ji„t = 160) of figure 2. 



With a negative interlayer coupling, the antiferromagnetic sublattice spins rotate 
with increasing temperatiu^e into the negative 2;-direction. For T > TpJ^gj, the 
magnetizations of both layers point in opposite directions, i.e. one has a ferrimagnetic 
situation. 

A remarkable behaviour can be seen in figure 2. With increasing interlayer 
couphng strength, the critical temperatures T^^^ — 26.24, 25.29, 22.00 and Tc = 
47.05, 46.34, 43.94 decrease in GFT. This is in contrast to the behaviour in MFT 
(Ti^,^^ = 51.19, 52.18, 55.86, Tc = 101.37, 103.28, 110.94), where the opposite is true 
(compare figures 2(a) and 2(b)). We attribute this surprising result to the quantum 
fiuctuations taken into account in Green's function theory but neglected in a mean field 
treatment. This behaviour depends, however, on the lattice type. It occurs for a sc 
lattice, whereas for fee stacking an increase of Tc with increasing interlayer coupling is 
obtained as in MFT. We show this by deriving a formula for the critical temperature 
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from GFT valid for both sc and fee bilayers. The magnetizations are eoUinearly arranged 
as T — > Tq, i.e. the x-components of the magnetizations vanish and the set of Green's 
functions reduces to two decoupled (2 x 2) problems from which one can obtain via 
the spectral theorem two equations for the magnetizations of both layers. Denoting the 
ratio of the magnetizations by ac = (•S'l)/ {SI) these equations become in this limit: 

^^•^TV ^ aia2 - «c(7int(k) Jint)^ " 2 ' ^^^^ 
2«c7c- E a^a2-ac(7int(k)Jint)2 ^ 2 ' ^^^^ 

with 

01 = go(-^i + Dl) + ac^int-'int - 7o(k) Ji, 

02 = acqo{J2 + Dl) + ^int Jint " Q;c7o(k) J2 ■ (28) 

For the notation see equations (8-12). The k summation is performed as an integral 
over the first Brillouin zone, and Tq is determined from a selfconsistent solution of the 
integral equations (26,27). Note that Tc is symmetric with respect to the sign of Jmt, if 
simultaneously ac —>■ —otc- 

In figure 4 the results for the critical temperatures for sc and fee stacking as a 
function of the interlayer coupling are shown. In the upper panel we see that for a sc 
lattice Tc decreases with increasing interaction strength in GFT, whereas it increases 
with increasing interlayer interaction strength in MFT. For fee stacking, Tc increases 
with increasing interlayer interaction strength for both GFT and MFT. This behaviour 
we attribute to a different action of the quantum fluctuations in sc and fee stacking. 
For large interlayer couplings the critical temperatures approach a saturation value in 
GFT, whereas they increase infinitely in MFT, which cannot be the correct behaviour. 

We now investigate the results of a sign change of the interlayer coupling Jint for 
the case of an sc lattice. In figure 5 we display the sublattice magnetizations for the 
ferromagnetic- antiferromagnetic bilayer at temperature T = as functions of the 
interlayer coupling Jinf An asymmetric behaviour of the magnetizations is observed 
with respect to the sign of Jinf The z-components of the sublattice magnetizations 
are positive for positive interlayer couplings and negative for negative couplings, while 
the x-components of the magnetizations of the ferromagnetic sublattices interchange 
their roles. More interesting is that the magnitudes of the corresponding components 
of the sublattice magnetizations are not identical for ±Jinf This situation with respect 
to the sign of the interlayer coupling does not change at finite temperatures. We 
discuss this without showing a figure. When the x-components of the magnetizations 
vanish above the temperature Tj^^^j, the magnetizations of the sublattices are collinear; 
for a positive interlayer coupling all z-components point in the positive 2;-direction, 
whereas for negative interlayer coupling the magnetizations of the antiferromagnetic 
sublattices point in the negative z-direction, i.e. opposite to the magnetizations of the 
ferromagnetic sublattices. The magnitudes of the magnetizations are larger for positive 
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Figure 4. Critical temperatures Tc{Ji-at) as a function of the interlayer coupling 
Jint of a simple cubic (upper panel) and a fee (lower panel) FM-AFM bilayer (Tc is 
normalised to rc(Jint = 0)). Results with mean field theory (MFT) and with Green's 
function theory (GET) are compared. Note the qualitatively different behaviour of the 
GFT results for sc and fee stacking (see the text). 



interlayer coupling. Even though the magnetization components are somewhat different 
for positive and negative interlayer couplings, numerical calculation and inspection 
of equations (26,27) show that the corresponding critical temperatures are identical. 
We attribute the observed asymmetries in the magnetizations in GFT to quantum 
fluctuations in the non-coUinear state, because there is no sensitivity on the sign of 
the interlayer coupling in MFT [6]. 

Finally, we discuss the influence of an external field on the sublattice 
magnetizations. As the only example we show in figure 6 the sublattice magnetizations 
at temperature T = as a function of a field in the direction of the AFM easy 
axis and a small perpendicular field component 5^ = 0.01. At = 0, one has a 
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GFT:FM-AFM Jfm=1 00 Jafm=-50 T=0 




W -0.5 



-200 -100 100 200 
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Figure 5. Sublattice magnetizations for the ferro- antiferromagnetic bilayer as 
function of the interlayer coupUng strength Ji„f 



GFT: FM-AFM bi-layer in external field 

JPM=1 00 J,p„=50 Ji„,=30 B^=0.01 T=0 




Figure 6. Sublattice magnetizations for the ferro- antiferromagnetic bilayer as funtion 
of a field in x-direction and a small field in z-direction. 

nearly orthogonal configuration of tlie FM (^-direction) and AFM (x-direction) layers. 
Owing to the positive interlayer coupling, both of the antiferromagnetic sublattice 
magnetizations have small positive z-components, and the ferromagnetic sublattice 
magnetizations have a small positive and a negative x-component. With increasing 
, a field-induced magnetic reorientation occurs, i. e. both of the x-components of 
the magnetization of the ferromagnet become positive and increase, whereas the z- 
components decrease. For very large fields, the magnetization of the ferromagnet points 
almost in the x-direction. Accordingly, for a small i?^-field the AFM layer preserves its 
nearly antiferromagnetic configuration and the almost orthogonal magnetic arrangement 
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with respect to the ferromagnetic layer. With increasing B^, the AFM sublattice 
magnetizations rotate more and more into the z-direction, with two small positive x- 
components. The situation changes only slightly if one puts the ^-component of the field 
to zero, B''^ = 0. In this case the main difference is that the reorientation transition is 
sharp at a reorientation field B^ = 3.79. The behaviour of the field dependence at finite 
temperatures is analogous, the only difference being that the magnetization components 
are reduced. 

4. Concluding remarks 

In the present paper we have developed a theory for coupling ferromagnetic and 
antiferromagnetic Heisenberg layers in the framework of many-body Green's function 
theory. The new feature is that we allow a non-coUinear orientation of the sublattice 
magnetizations for both of the ferromagnetic and antiferromagnetic layers. We apphed 
the theory to bilayer systems with an in-plane orientation of the magnetizations. 
The theory yields results which are in many instances qualitatively similar to 
results previously obtained by mean field theory [6]. Owing to missing quantum 
fluctuations, MFT cannot, however, give the suppression of the magnetizations of the 
antiferromagnetic sublattices at low temperature. The main difference between MFT 
and GFT is, as for the uncoupled systems, the very different temperature scale. 

An interesting effect is observed for the critical temperature Tq for systems with 
JpM > |<^afm|- For a sc lattice, Tq decreases with increasing interlayer coupling J^^i 
in GFT, contrary to the behaviour in MFT, where the critical temperature increases 
with increasing Jint. For fee stacking, however, for both GFT and MFT the critical 
temperature increases with increasing Jint- We attribute the different behaviour of the 
GFT results to a different action of the quantum fluctuations for sc and fee stacking. 

In the future we will study multilayer systems with mixed FM and AFM couplings. 
The code is written in such a way that it also allows the description of the coupling of 
ferromagnetic layers by an antiferromagnetic interlayer coupling and vice versa. Because 
we have already included magnetic fields, the theory may also be the basis for studying 
the exchange bias effect, where it seems, however, to be necessary to include interface 
disorder [21] in some way, for instance by introducing more sublattices per layer with 
different magnetic arrangements. 
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